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Abstract  The  far  infrared  spectra  of  (100),  (010),  and 
(OOl)-oriented  RDX  single  crystals  were  measured  as  the 
crystal  was  rotated  about  the  axis  perpendicular  to  the 
polarization  plane  of  the  incident  radiation.  Absorption 
measurements  were  taken  at  temperatures  of  both  20  K  and 
295  K  for  all  rotations  using  terahertz  time-domain 
spectroscopy.  A  number  of  discrete  absorptions  were  found 
ranging  from  10-100  cm-1  (0.3-3  THz).  The  absorptions 
are  highly  dependent  on  the  orientation  of  the  terahertz 
polarization  with  respect  to  crystallographic  axes. 

Keywords  IR  spectroscopy/Raman  spectroscopy  • 

Organic  compounds/trace  organic  compounds  • 
Spectroscopy/instrumentation 

Introduction 

Terahertz  time-domain  spectroscopy  (THz-TDS)  shows 
potential  as  a  noninvasive  technique  for  the  detection  of 
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explosives.  Common  energetic  materials  such  as  RDX, 
HMX,  PETN,  and  TNT  have  unique  spectral  signatures  in 
the  <100  cm-1  range  [1-15].  At  visible  wavelengths,  many 
explosives  look  remarkably  similar  to  common  baking 
ingredients  such  as  sugar  or  modeling  clay  but  have  unique 
signatures  in  the  THz  range.  Light  in  the  THz  region  (e.g., 
below  100  cm-1)  will  pass  through  visibly  opaque  materials 
such  as  clothing,  luggage,  paper,  etc.  Thus,  THz  radiation 
can  be  used  to  positively  identify  unknown  substances 
insides  opaque  materials — a  highly  useful  trait  for  security 
applications  [3,  16,  17]. 

While  well-characterized  [1,  6,  9],  the  spectrum  of  RDX 
below  100  cnT1  is  poorly  understood.  Measurements  on 
powdered  RDX  give  an  absorption  spectrum  that  is  an 
average  over  many  randomly  oriented  crystals.  To  compli¬ 
cate  matters,  the  manufacturing  process  has  some  effect  on 
the  convoluted  spectrum  as  well.  Absorption  spectra 
measured  on  RDX  from  different  manufacturers  show 
significant  variation  [18].  Assignment  of  absorptions  to 
the  vibrational  modes  of  the  crystal  structure  is  difficult 
when  using  powders  because  the  high  density  of  absorp¬ 
tions  make  distinguishing  each  individual  peak  difficult  and 
each  absorption  cannot  be  assigned  to  a  given  crystal 
orientation.  Calculations  on  isolated  RDX  molecules  are 
not  able  to  properly  account  for  the  interactions  measured 
in  the  solid  phase  as  the  majority  of  the  absorptions  seen  in 
the  solid  phase  are  intermolecular  in  character  [19]. 
Progress  has  been  made  recently  into  the  calculated  spectra 
of  HMX  and  PETN  in  the  solid  phase  [20-22],  but  RDX  is 
a  more  difficult  computational  problem  due  to  the  large 
number  of  molecules  in  the  unit  cell  (Z=  8  for  oc-RDX  [23], 
Z=2  for  (3 -HMX  [24],  Z=2  for  PETN  [25]).  Recent  work 
by  Allis  et  al.  [26]  has  produced  reasonable  agreement  with 
powdered  RDX. 
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Developing  models  that  can  correctly  characterize  the 
complex  absorption  patterns  of  RDX  in  the  <100  cm-1  region 
is  critically  important.  This  phonon-rich  region  of  the 
spectmm  contains  information  relevant  to  the  validation  of 
equation  of  state  parameters  used  to  model  the  detonation 
process.  Accurate  modeling  would  help  clarify  the  role  these 
phonon  modes  play  in  shock-initiation  of  explosives.  It  has 
been  postulated  that  energy  is  tunneled  from  the  shock  into 
individual  molecules  through  multiphonon  up-pumping  via 
delocalized  low-frequency  (<6  THz)  phonons  [27-30],  but 
little  experimental  data  on  these  phonon  modes  exist. 

In  an  effort  to  resolve  and  identify  the  vibrational  modes 
of  RDX,  single  crystals  of  a-RDX  have  been  studied  via 
THz-TDS.  This  work  is  a  continuation  of  the  study 
published  by  Barber  et  al.  in  2005  [19].  Crystals  cut  at 
three  orthogonal  orientations  were  mounted  and  oriented 
via  X-ray  diffraction.  An  azimuthal  rotation  stage  with  a 
center  hole  for  performing  transmission  spectroscopy  was 
mounted  to  a  cryostat,  enabling  spectra  to  be  collected  at 
~20  K.  The  results  of  this  study  show  terahertz  absorptions 
not  previously  observed  from  powders  or  formulations. 

Experimental  section 

For  these  measurements,  five  different  RDX  samples  were 
used — two  (lOO)-oriented  samples,  two  (OlO)-oriented 
samples,  and  one  (OOl)-oriented  sample.  The  samples  were 
cut  from  large  single  crystals,  nominally  5  cm  on  a  side. 
UK-manufactured  RDX  was  dried,  purified  by  Soxhlet 
extraction  in  acetone  and  recrystallized  in  acetone  prior  to 
use.  The  crystals  were  grown  by  slow  evaporation  of 
acetone  at  32  °C,  in  a  1 ,000  mL,  spoutless  beaker  covered 
by  a  Parafilm-coated  watch  glass.  The  bottom  third  of  the 
beaker  was  immersed  in  a  heating  bath  and  the  top  was 
open  to  air  at  room  temperature.  This  method  allowed  for 
convective  mixing  and  the  reflux  continuously  washed  the 
beaker  sides  preventing  parasitic  growth.  Seeds  were  grown 
from  acetone  by  a  similar  procedure. 

Contact  goniometry  was  used  to  locate  planes  in  RDX 
by  redundant  reference  to  crystal  facets.  The  samples  were 
cut  using  a  low-speed  diamond  impregnated  wire  saw  using 
a  solution  of  Alconox  in  deionized  water  as  a  lubricant. 
Samples  were  measured  in  an  as-cut  configuration;  no 
further  polishing  was  performed.  RDX  crystals  are  ex¬ 
tremely  fragile,  limiting  the  minimum  sample  thickness  that 
can  be  cut  without  fracture  to  approximately  1  mm.  The 
(001)-  and  (OlO)-oriented  samples  were  nominally  1  mm 
thick,  the  (100)  sample  was  1.5  mm  thick.  All  samples  were 
approximately  1  cmxl  cm  wide.  The  orientations  of  the 
crystallographic  axis  and  the  orientation  of  the  in-plane 
axes  were  verified  by  transmission  Laue  X-ray  diffraction. 
Typically,  the  uncertainty  in  axis  orientation  is  within  one- 


half  degree  of  the  reported  face.  Plane  indices  are  given 
with  respect  to  the  Pbca  space  group  of  RDX  [23].  The 
structure  of  the  RDX  crystal  cut  along  the  (100),  (010),  and 
(001)  faces  is  shown  in  Fig.  1. 

The  THz-TDS  system  was  similar  to  those  used  in  previous 
work  [31]  with  some  modifications  to  accommodate  the 
specialized  helium  cryostat.  Photoconductive  antennas 
backed  with  silicon  lenses  served  to  both  generate  and  detect 
the  collimated  THz  radiation,  while  confocal  optics  guided 
the  THz  beam  through  the  samples.  The  beam  was 
achromatically  focused  to  a  spot  about  6  mm  in  diameter  at 
the  sample  face.  The  maximum  spectral  resolution  of  the 
system  was  1.1  cm-1.  The  polarization  purity  was  nominally 
100:1,  which  is  typical  for  photoconductive  switch  systems. 

The  cryostat  and  sample  mounting  structures  were  specially 
made  to  allow  sample  rotation  and  cooling  under  vacuum. 
Samples  were  mounted  on  the  front  end  of  the  rotating  stage, 
essentially  a  copper  tube  through  which  the  focused  THz  beam 
propagated.  The  tube  diameter  was  sufficiently  large  to  avoid 
interactions  with  the  THz  beam.  A  ring  gear  was  attached  to 
the  back  of  the  tube  and  engaged  to  a  pinion  gear  that  was 
magnetically  coupled  through  the  vacuum  chamber  to  a 
computer-controlled  stepper  motor.  This  allowed  accurate 
rotational  stepping  of  the  sample  orientation  under  vacuum. 
The  tube  was  thermally  connected  to  the  main  cryostat  cold 
finger  by  a  flexible  copper  band.  The  thermal  link  was 
sufficient  to  reduce  the  sample  temperature  to  18  K  without 
interfering  with  sample  rotation.  Sample  temperatures  were 
monitored  by  a  silicon  sensor  attached  adjacent  to  the  sample. 

Initial  mounting  of  the  samples  into  the  cryostat  and 
orientation  of  the  rotation  stage  introduced  an  angular 
inaccuracy  within  approximately  +/-3°.  Samples  were 
cooled  in  a  helium  cryostat  at  0.025  K/s  to  avoid  cracking. 
Even  at  these  rates,  one  each  of  the  (010)  and  (100)  samples 
cracked.  Comparisons  of  the  absorptions  of  the  cracked  and 
uncracked  samples  show  that  the  cracks  produced  no 
significant  variations  in  the  data.  Over  the  full  1 80°  rotation, 
the  sample  temperature  ranged  from  18-23  K  and  had  a 
maximum  variation  of  0.5  K  during  any  one  scan.  The 
variation  in  temperature  was  caused  by  the  varying  thermal 
link  of  the  rotation  stage  as  it  changed  position.  Based  on  the 
temperature-dependent  absorption  measurements  of  RDX 
by  Melinger  [15],  Barber  [19],  and  Burnett  [32],  the  shape 
and  position  of  the  spectra  are  expected  to  remain 
unchanged  over  this  temperature  variation. 

Analysis 

The  time-domain  data  were  acquired  for  each  sample 
orientation  over  a  180°  rotation  around  the  primary 
crystallographic  axis  using  angle  increments  of  9°  at  22  K 
and  22.5°  at  295  K.  Reference  data  sets  with  the  sample 
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Fig.  1  Representation  of  RDX  along  the  primary  directions  a  100  b  010  c  001.  Carbon  is  shown  in  black ,  nitrogen  in  blue ,  hydrogen  in  green , 
and  oxygen  in  red 


moved  out  of  the  path  of  the  beam  were  taken  at  the 
beginning  and  end  of  each  sample  data  set.  The  ratio  of  the 
Fourier-transformed  THz-TDS  sample  data  to  a  reference 
data  set  yielded  the  complex  transmission  function.  Figure  2 
shows  the  results  of  a  typical  terahertz  reference  amplitude. 
Also  shown  are  the  transmission  amplitude  through  RDX 
(010)  with  the  polarization  oriented  along  the  oaxis  and  the 
noise  level  of  the  detector/lock-in  amplifier.  The  system  is 
capable  of  producing  measurable  quantities  of  radiation  to 
150  cm-1.  The  combination  of  decreased  emission  ampli¬ 
tude  above  -100  cm-1  and  strong  sample  absorption 
resulted  in  noisy,  inconsistent  data.  None  of  the  data  above 
100  cm-1  is  discussed  in  this  paper,  but  the  data  to 
135  cm-1  is  included  in  the  data  files  in  the  electronic 
supplemental  material.  Strong  absorption  features,  such  as 
the  57  cm-1  feature  shown  in  Fig.  2,  are  completely 
absorbing.  Since  the  minimum  transmission  is  not  measur¬ 
able,  the  data  produces  an  artificial  maximum  value  of  the 
absorption  index  k  or  the  power  absorption  coefficient  a. 

Since  the  samples  are  homogenous  single-crystal  specimen 
with  planar  faces  and  known  thicknesses,  indices  of  refraction 
and  absorption  values  as  a  function  of  crystal  orientation  can 
be  extracted.  Initially,  a  Beer’s  law  analysis  was  used  to 
determine  measured  absorption  coefficients  (not  correcting  for 
Fresnel  reflections)  through  the  sample  at  each  orientation. 
This  was  refined  by  extracting  the  complex  index  of  refraction 
n  and  absorption  index  n  from  the  THz-TDS  data  using  the 
method  of  Duvillaret  et  al.  [33].  The  similarity  of  the  two 
methods  provided  a  check  that  the  fitting  routines  were 
working  correctly.  The  fits  were  also  verified  by  mnning  the 
extracted  results  through  the  forward  propagation  equation, 
which  matched  the  measured  data. 


Results  and  discussions 

Figures  3-4  show  the  angular  dependence  of  k-  and  n- 
values  for  samples  of  RDX  (100),  (010),  and  (001)  taken  at 


20  K.  The  data  are  presented  as  radial  contour  plots  with 
the  THz  absorption  spectrum  extending  radially  outward 
from  the  center  to  a  maximum  of  100  cm-1.  Dashed  circles 
indicate  20  cm-1  increments.  The  values  of  n-  or  k-  are 
plotted  as  a  function  of  rotation  about  the  sample.  Principle 
in-plane  axes  are  indicated  on  the  plots.  Since  data  were 
taken  over  180°  of  rotation  instead  of  a  full  360°  rotation, 
the  data  were  reproduced  over  the  180-360°  range  to 
generate  a  full  rotation  spectrum.  RDX  has  C2  symmetry  so 
data  from  0-180°  should  be  identical  to  180-360°.  The 
reader  should  take  some  caution  with  regions  of  the 
spectrum  showing  strong  absorption.  These  absorptions 
are  strong  enough  to  completely  block  all  incident  THz 
radiation,  so  that  the  absorption  values  in  these  locations 
should  be  considered  as  a  minimum  instead  of  exact. 

All  three  orientations  of  RDX  in  Fig.  3  show  a  complex, 
feature-rich  absorption  structure.  Significant  absorption 
features  extend  from  as  low  as  12  cnT1  to  well  above 
100  cm-1.  Below  -80  cm-1,  RDX  has  a  number  of 


Fig.  2  Electric  field  amplitude  of  the  THz  spectral  emission  with  no 
sample  {black),  through  (010)  RDX  with  the  polarization  oriented 
along  the  c-axis  {red),  and  the  noise  level  of  the  detector/lock-in 
amplifier  {green) 
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Fig.  3  Comparison  of  the  absorption  index  k  of  RDX(100),  RDX 
(010),  and  RDX(001)  over  the  7-100  cm-1  range  at  T=20  K.  The 
dashed  circles  correspond  to  20  cm-1  wavenumber  intervals. 
Crystallographic  axes  are  shown  as  the  labeled  solid  lines 


rotationally  dependent  absorptions  strongly  correlated  to  its 
crystallographic  orientation  with  respect  to  the  THz 
polarization,  above  ~80  cm"1,  absorptions  are  generally 
rotationally  invariant. 

Broad  absorption  features  below  80  cm-1  are  found 
along  the  primary  crystallographic  axes.  The  Z?-axis  has  a 
broad  absorption  from  27-42  cm-1  in  (OOl)-oriented  RDX 
and  from  48-60  cm-1  in  (lOO)-oriented  RDX.  The  a- axis 
has  a  broad  absorption  from  55  cm-1  to  59  cm-1  seen  in 
both  the  RDX(010)  and  RDX(001)  orientations.  The  oaxis 
has  a  broad  absorption  from  25-40  cm-1  in  the  (100) 
orientation,  although  there  is  no  corresponding  oaxis 
absorption  in  the  (OlO)-oriented  material.  In  addition,  there 
are  two  broad  low-frequency  absorptions  located  at  45° 
between  in-plane  axes.  RDX(100)  shows  an  absorption  at 
18  cm-1  located  between  the  b-  and  c-  axes  and  RDX(010) 


shows  an  absorption  at  12  cm-1  located  between  the  a-  and 
Z?-axes. 

Narrow,  rotationally  invariant  absorptions  with  1  cm-1 
FWHM  are  found  at  40  cm-1  in  RDX(100),  55  cm-1  in 
RDX(010),  and  at  25  cm-1  in  RDX(001).  Narrow, 
orientation-dependent  absorptions  associated  with  the  o 
axis  are  seen  in  the  (100)  and  (010)  orientations  at  27  cm-1. 
Other  narrow  absorptions  are  roughly  located  equally 
between  the  primary  in-plane  axes  at  24  cm-1,  48  cm-1, 
67  cm-1  in  RDX  (001)  and  40  cm-1  and  45  cm-1  in  RDX 
(010).  A  summary  of  the  in-plane  absorptions  found  in 
RDX  is  listed  in  Table  1 . 

The  absorbing  states  below  40  cm-1  in  RDX  are  unique 
when  compared  to  HMX  [34]  or  PETN  [34].  HMX  and 
PETN  do  not  show  absorptions  at  such  low  frequencies. 
The  lowest  absorption  for  PETN  is  55  cm-1  and  for  HMX 
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Fig.  4  Comparison  of  the  index  of  refraction  of  RDX(100),  RDX 
(010),  and  RDX  (001)  from  7-100  cm-1  at  T=20  K.  Crystallographic 
axes  are  shown  as  the  labeled  solid  lines 
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Table  1  Summary  of  absorption  found 

in  (100)-,  (010)-,  and  (OOl)-oriented  RDX 

Orientation 

Frequency  cm'1  22K 

Frequency  cm'1  295K 

Absorption  Index 

oo 

Absorption  coefficient 
oc  (cm-1) 

Relation  to  axes 

(100) 

17.6 

14.2 

0.19 

42 

45°  between  6,c 

25.4 

0.14 

45 

c- axis 

27-37 

22-37 

0.23 

78 

c-axis 

39.5 

0.07 

35 

R.I 

46-61 

40-58 

0.15 

87-115 

6-axis 

66.7-77.3 

>60 

0.12 

>75 

R.I 

83.3-87.5 

0.09 

>94 

R.I 

94-97 

0.09 

>106 

R.I 

(010) 

25 

19.8 

0.05 

16 

c-axis 

38.3 

38.8 

0.10 

48 

30°  from  a-axis 

42.2 

0.10 

53 

30°  from  a-axis 

46-60 

42-53 

0.18 

104-135 

c-axis 

51.54 

48.4 

0.15 

91-101 

R.I 

64.3 

0.10 

81 

a-axis 

67-77 

66-72 

0.13 

>110 

a-axis 

75.5 

0.06 

>56 

R.I 

82-86 

0.08 

>82 

R.I 

>94 

0.08 

>94 

R.I 

(001) 

13.8 

10.3 

0.18 

31 

45°  between  a,b 

23.2 

0.08 

23 

30°  from  6-axis 

26.5 

18.9 

0.05 

17 

R.I.  and  30°  from  6-axis 

28.41 

20-38 

0.14 

50-72 

6-axis 

48.3 

43.5 

0.04 

24 

45°  between  a,  b 

51-56 

46-51 

0.09 

58-63 

a-axis 

65.3 

58-67 

0.03 

>46 

45°  between  a ,  b 

68.9 

0.05 

>43 

15°  from  6 -axis 

69-75 

0.07 

>61 

6-axis 

75-79 

>68 

0.07 

>66 

R.I. 

82-89 

0.07 

>72 

R.I. 

93.2 

0.04 

>47 

30°  from  6-axis 

94-99 

0.04 

>47 

6-axis 

>97 

0.04 

>49 

R.I. 

Rotationally  invariant  absorptions  are  listed  as  ‘R.I.’  in  the  table 


is  40  cm-1.  The  causes  of  the  low-frequency  absorption  in 
RDX  are  thought  to  be  related  to  the  six  extra  molecules  in 
the  unit  cell  compared  to  PETN  and  HMX  and  the  larger 
unit  cell  of  RDX  (13.192-11.574-10.709  A  for  RDX  [23] 
compared  to  9.380-9.380-6.71  A  for  PETN(I)  [25]  and 
6.54-11.05-8.7  A  for  HMX  [24]).  The  absorption  peaks  of 
RDX  in  the  25-30  cm-1  range  are  particularly  significant 
because  they  are  responsible  for  the  ~27  cm-1  “fingerprint” 
peak  of  powdered  RDX.  The  absorptions  of  single-crystal 
RDX  below  27  cm-1,  namely  the  broad  peaks  located  at 
18  cm-1  and  12  cm-1  in  RDX(100)  and  (001),  respectively, 
have  not  been  previously  reported.  In  the  powdered 


spectrum,  the  27  cm-1  absorption  is  a  broad  dominant 
feature  that  masks  the  12  cm-1  and  18  cm-1  absorptions. 

Figure  4  shows  the  frequency-dependent  index  of 
refraction  of  single-crystal  RDX  as  a  function  of  orienta¬ 
tion.  The  index  of  refraction  increases  from  ~1.6  to  ~3  over 
the  frequency  range  of  7-100  cm-1.  RDX,  known  to  be 
birefringent  in  the  visible  region  [35],  is  also  birefringent  in 
the  7-100  cm-1  range.  Symmetry  points  in  the  angular 
variation  of  the  index  of  refraction  correspond  identically 
with  the  primary  in-plane  axes. 

This  birefringence  can  be  used  to  find  the  principle 
in-plane  axes  of  RDX  using  THz-TDS.  The  time  delay 
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in  the  onset  of  the  THz-TDS  data  varies  as  a  function  of 
orientation.  The  extremes  of  the  measured  time  delay 
give  the  two  in-plane  axes.  The  minimum  time  delay 
corresponds  to  the  axes  with  the  smaller  indices  of 
refraction  and  the  maximum  time  delay  corresponds  to 
the  larger  indices  of  refraction.  For  example,  the 
minimum  time  delay  in  the  pulse  envelope  of  RDX 
(100)  would  locate  the  Z?-axis  and  the  maximum  delay 
would  determine  the  oaxis. 

Figure  5  shows  a  comparison  of  the  power  absorption 
coefficients  of  the  three  orthogonal  orientations  taken  at 
295  K  and  20  K.  The  absorption  coefficients  were 
determined  through  Beer’s  law  analysis  of  the  data  and 
have  not  been  corrected  for  Fresnel  reflections.  The 


absorption  spectra  of  RDX  at  295  K  are  very  similar  to 
the  absorption  spectra  at  22  K  for  all  orientations.  Upon 
cooling  the  samples,  the  spectral  features  become  sharper 
and  shift  to  higher  frequencies.  No  other  major  changes  in 
the  absorption  are  found.  A  summary  of  the  frequency 
shifts  due  to  temperature  are  included  in  Table  1. 

This  is  the  first  work  to  present  the  rotational  dependence 
of  the  absorptions  of  RDX  single  crystals,  so  comparison 
with  other  published  work  is  not  possible.  We  measured 
two  different  RDX(OIO)  and  RDX(IOO)  samples,  and 
performed  two  different  measurements  on  the  same  RDX 
(001)  sample  to  ensure  that  the  data  were  reproducible. 
There  were  no  major  sample-to-sample  or  measurement-to- 
measurement  differences  in  the  absorption  spectra.  This 


Fig.  5  Comparison  of  the  pow¬ 
er  absorption  coefficient  oc  for 
RDX(100)  at  295  K  (left)  and  at 
20  K  (right)  for  a  RDX(100),  b 
RDX(010),  and  c  RDX(001) 
over  the  7-100  cm-1  spectral 
range.  The  dashed  circles  corre¬ 
spond  to  20  cm-1  wavenumber 
intervals.  The  contour  intervals 
give  the  absorption  coefficient 
in  units  of  cm-1 


■  ooo 

■  9.00 

■  J  3.00 

■  27.00 

■  36.00 

■  45.00 

■  MOO 
63.00 

■  72.00 

■  81.00 
■  90-00 


■  0.00 
■  8.00 
■  16.00 

■  24.00 

■  32.00 

■  40.00 

■  48.00 

■  56.00 

■  64.00 

■  72.00 

■  80.00 


■  0,00 
■  900 
13.00 

■  2700 

■  36.00 
45.0O 
54  00 

■  63.00 

■  72.00 

■  81,00 
■  90.00 
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study  does  overlap  with  the  work  of  Barber  et  al. 
presenting  the  absorption  spectrum  of  RDX(IOO)  along 
a  single  random  orientation  from  20-295  K  [19]. 
Comparing  that  spectrum  with  the  RDX(IOO)  absorptions 
shown  above  indicates  that  the  sample  used  by  Barber  et 
al.  was  roughly  oriented  along  the  Z?-axis  in  relation  to  the 
THz  polarization. 

An  averaged  absorption  spectrum  of  all  the  angular 
contributions  can  be  compared  to  powder  spectra  if  we 
assume  that  the  powder  spectrum  is  composed  of  an  equally 
weighted  contribution  from  all  possible  crystal  orientations. 
To  generate  the  equivalent  of  a  powder  absorption 
coefficient,  we  first  calculate  the  average  measured  light 
intensity  as  a  function  of  the  various  angles  measured  from 
all  three  crystal  orientations: 

n 

,T\  J2ExP(~aet) 

/  1  \  e  /!\ 


where  oie  is  the  absorption  coefficient  from  a  given  angular 
contribution,  t  is  the  thickness  of  the  sample,  and  n  is  the 
number  of  angular  contributions  measured.  We  are  able  to 
calculate  this  averaged  intensity  using  the  values  of  k  from 
Fig.  3,  and  relating  the  absorption  coefficient  a  to  the 
absorption  index  k  and  wavelength  A  by:  a=4nK/\.  Once 
we  have  the  average  measured  light  intensity  from  an  equal 
distribution  of  angles,  then  we  can  calculate  an  average 
absorption  coefficient  by  solving  Beer’s  law: 


Figure  6  shows  the  results  of  this  analysis  using  the 
values  of  k  determined  from  all  three  orientations  of  RDX 
at  22  K.  Care  should  be  taken  when  viewing  the  averaged 


Fig.  6  Average  spectrum  of  single  crystal  RDX  at  20  K  composed  of 
an  equally  weighted  average  of  all  the  angular  contributions 


spectrum  because  the  values  of  absorption  index  k 
determined  for  the  strongest  absorptions  in  Fig.  3  are 
artificially  truncated  due  to  dynamic  range  limitations,  so 
their  contributions  to  the  average  spectrum  will  be 
artificially  lowered.  Nevertheless,  strong  absorptions  are 
found  at  27,  40,  51-53,  75,  83,  and  96  cm-1.  Allis  et  al. 
[26]  and  Burnett  et  al.  [32]  measured  the  spectrum  of 
powdered  RDX  at  5  K.  There  are  similarities  between  our 
calculated  powdered  spectrum  and  the  measured  spectrum; 
however,  they  are  not  identical.  Perhaps  the  measured 
powder  spectrum  is  not  a  true  representation  of  all  crystal 
orientation.  Crystal  morphology  may  cause  the  crystals  to 
pack  along  given  orientations.  We  also  do  not  know  how 
the  surface  effects  contribute  to  the  total  RDX  powder 
spectrum.  Surface-related  spectral  features  will  be  much 
more  prominent  when  using  powders  compared  to  single 
crystals. 

Melinger  et  al.  [15]  measured  high  resolution  absorp¬ 
tion  of  ordered  poly  crystalline  RDX  films  at  1 3  K,  finding 
19  separate  peaks  between  0-116  cm-1.  It  is  unclear  how 
those  peaks  compare  to  this  work  or  to  the  powdered 
spectrum  of  Burnett.  If  the  films  consisted  of  preferential¬ 
ly  oriented  polycrystalline  RDX,  as  the  authors  suggest, 
we  could  not  decipher  a  single  orientation  of  the  RDX 
films  necessary  to  compare  with  the  data  in  this  paper. 
Perhaps  the  thin  film  sample  was  composed  of  several 
orientations  with  respect  to  the  surface  and  the  THz 
polarization. 

Conclusions 

We  present  absorption  spectra  of  RDX  (100),  (010)  and 
(001)  from  7-100  cm-1  at  295  K  and  20  K  as  a  function 
of  orientation  with  respect  to  THz  polarization.  Results 
demonstrate  the  expected  dependence  of  spectral  features 
on  angular  orientation  through  alignment  of  the  oscillat¬ 
ing  dipole  to  the  THz  field.  Absorptions  as  low  as 
12  cm-1  have  been  observed,  the  lowest  frequency 
absorptions  reported  for  RDX  to  date.  Line  widths  are 
predominantly  broad  (10  cm-1  FWHM),  with  some  modes 
as  narrow  as  1  cm-1.  It  should  be  possible  to  make 
assignments  for  the  spectrally  resolved  absorptions  by 
comparing  these  results  with  normal  mode  calculations  of 
an  a-RDX  crystal. 
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